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We study the symmetry classes of graphene quantum dots, both open and closed, through the 
conductance and energy level statistics. For abrupt termination of the lattice, these properties 
are well described by the standard orthogonal and unitary ensembles. However, for smooth mass 
confinement, special time-reversal symmetries associated with the sublattice and valley degrees of 
freedom are critical: they lead to block diagonal Hamiltonians and scattering matrices with blocks 
belonging to the unitary symmetry class even at zero magnetic field. While the effect of this structure 
is clearly seen in the conductance of open dots, it is suppressed in the spectral statistics of closed 
dots, because the intervalley scattering time is shorter than the time required to resolve a level 
spacing in the closed systems but longer than the escape time of the open systems. 
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Single atomic layers of graphite, known as graphene, 
have attracted intense experimental and theoretical at- 
tention due to its unusual band structure and hence ex- 
otic electronic properties P, Moreover, graphene's 
true two-dimensional nature and high mobility make it 
an attractive alternative for studying low dimensional 
electron systems such as quantum dots d, 0, [^, @1 ■ Re- 
cent experiments on the spectra of graphene quantum 
dots 3 found evidence for a time reversal (TR) symme- 
try broken state in the absence of magnetic field, raising 
questions about the possible origin of such states. Some 
time ago Berry and Mondragon [3] proposed one such 
mechanism of TR symmetry breaking, namely infinite- 
mass confinement. In graphene dots, edge magnetization 
might produce such an effective mass term at the edges of 
the graphene flakes [1, Q , but whether this term is strong 
enough to change the universality class of the graphene 
quantum dots has not been established. 

In this work we study universalities in the spectrum 
and conductance of graphene quantum dots in both the 
closed Coulomb blockade and the open ballistic regime, 
respectively. Universal properties are generally deter- 
mined by the symmetries of the Hamiltonian or the scat- 
tering matrix [ld | . Thus one expects that the universality 
class displayed by the spectrum of a closed quantum dot 
should be identical to that displayed by the conductance 
of a corresponding open dot. Here we show that this 
naive expectation is not true: the universality class of the 
conductance can be different from that of the spectrum. 
The main reason behind this paradox is the separation of 
time scales characterizing the conductance (escape time) 
and the spectrum (Heisenberg time, i.e. inverse level spac- 
ing), allowing scattering times to be smaller than one but 
larger than the other. To demonstrate this scenario, we 
first focus on closed graphene dots and show that their 
spectral statistics is described by the orthogonal symme- 




FIG. f: Systems studied numerically (schematic). (a),(b) 
Africa billiard. (c),(d) Half-stadium with two identical leads; 
left-right symmetry is broken by cutting out circular segments 
at the top left and bottom right. The graphene lattice is ter- 
minated abruptly in (a) and (c), while smooth mass confine- 
ment is used in (b) and (d). 



try class even in the presence of collinear edge magne- 
tization, ruling out Berry and Mondragon's mechanism 
7] for TR symmetry breaking in this case. We next 
treat quantum transport through open dots and show 
that edge magnetism is enough to change the symmetry 
class, so that the conductance is described by the unitary 
ensemble. 

Symmetries of the Hamiltonian — The effective Hamil- 
tonian for low energies and long length scales is the well- 
known Dirac Hamiltonian (the spin is omitted here). 
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where the Pauli matrices ai and act on sublattice 
and valley degrees of freedom, respectively, and the in- 
dex i = denotes the unit matrix. The boundary of 
the graphene flake is critical for its properties; we dis- 
tinguish two physically relevant boundary types: (i) an 
abrupt termination of the graphene lattice, and (ii) con- 
finement by the mass term in Eq. (1). In the former case, 
m{x, y) = 0; the boundary is disordered on the lattice 
scale and contains valley mixing armchair edges. In case 
(ii), while the lattice eventually terminates, the confine- 
ment is due to the smooth mass term which prevents the 
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particles from feeling the rough boundary and thus sup- 
presses the intervalley scattering. The mass term may 
originate from an effective staggered potential caused by 
possible edge magnetization of graphene flakes d, . 

The symmetries of the problem are defined through 
three antiunitary operators 11, T^: time reversal T, and 



two "special time reversal" operators %i and %, associ- 
ated with either the sublattice or valley pseudospin: 

r=(cro®r^)C, -i{ay(^To)C, % = -i{ao®Ty)C. 

(2) 

C denotes complex conjugation. For abrupt termination, 
the two sublattices are inequivalent and boundary scat- 
tering mixes the valleys, so both special TR symmetries 
are irrelevant [31 . For B = 0, T commutes with HcS, 
leading to the orthogonal symmetry class. When B ^ 0, 
the Hamiltonian falls into the unitary ensemble. 

For smooth mass confinement, intervalley scattering is 
small, so that the system largely consists of two indepen- 
dent subsystems, one for each valley. Each subsystem 
lacks TR symmetry, even at zero magnetic field, because 
T commutes only with the full i?off, while %i is broken 
by the mass term. Thus, the Hamiltonian for a single val- 
ley corresponds to the unitary symmetry class. For zero 
magnetic field, however, commutes with 7^ while 
= —I. Kramers' theorem then guarantees the de- 
generacy of the eigenvalues of the full Hamiltonian 14 1 . 



Since the Ty in 7^ switches the valleys, the degenerate 
states do not lie in the same valley. Thus the Hamiltonian 
consists of two degenerate blocks with unitary symmetry. 
Upon applying a magnetic field, Hcs does not commute 
with and the valleys are no longer degenerate. 

Spectral statistics — To exhibit the universality classes 
of closed graphene dots, we focus on the level spacing 
distribution for an Africa billiard with either abrupt 
termination or smooth mass confinement (Fig.[T]). For the 
numerical work, we use the tight-binding Hamiltonian 
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where i and j are nearest neighbors. The staggered 
potential rrii = m{xi,yi), corresponding to a mass 
term, is positive (negative) if i belongs to sublat- 
tice A (B). A magnetic field can be introduced via 

tij — — i exp JT;' A • drj , with the flux quantum 

^0 = h/ e. The lattice points are determined by cutting 
an Africa billiard out of a graphene plane [Fig.[TJa),(b)] 
with X being a zigzag direction. For smooth mass con- 
finement, the mass term is zero in the interior but 
non-zero within a distance W of the boundary [see 
Fig.[ljb)]; it starts from zero at the inner border of this 
region (black line in sketch) and increases quadratically: 



,y) 



[5{x, y) — W] /2, where S{x, y) is the dis- 



tance to the boundary and w is a constant. 

Fig. [21 shows the level spacing distribution for both 
abrupt termination and smooth confinement in an Africa 




FIG. 2: (Color online) Level-spacing distribution P{S) for 
an Africa flake consisting of 68169 carbon atoms using about 
3000 energy levels in the range [—0.5 0.5 t\. (a),(b) Abrupt 
lattice termination. (c),(d) Smooth mass confinement with 
= 0.15 y/t/a. (e),(f) Smooth mass confine- 
ment with W = 16.5V3a, ui = 0.041 y/t/a. a « 0.25 nm is 
the graphene lattice constant. The left panels are for $ = 0, 
while the right panels are for $ = 0.7"l?o- Insets in each panel 
present the integrated distributions C{S) = P{S')dS' . 
Numerical results are shown with solid thick black lines, 
whereas the thin lines are for Poisson (green solid) , GOE (red 
dotted), and GUE (blue dashed) statistics. 



graphene dot. For abrupt termination (top panels in 
Fig-Hi 7 the statistics are consistent with the Gaussian 
orthogonal ensemble (GOE) when _B = and with the 
Gaussian unitary ensemble (GUE) upon introduction of 
a magnetic field. This is expected from the symmetry 
considerations above. 

For smooth mass confinement, the results are surpris- 
ing: the statistics are not the expected GUE but rather 
are GOE for large systems [Fig.[2Kc)] with a crossover 
to Poisson for smaller systems [Fig.[2Je)]. This crossover 
reflects the role of localized edge states present for ener- 
gies near the Dirac point which follow Poisson statistics. 
Edge states dominate in small systems, but for larger sys- 
tems their spectral weight diminishes, giving rise to the 
crossover to GOE statistics. We believe this is why the 
numerical level statistics in Rcf. fl^ does not fit well to 
either Poisson or Gaussian ensembles. 

The reason that we find orthogonal rather than uni- 
tary statistics for a large dot is more subtle: Though our 
mass confinement is fairly smooth, there is some residual 
intervalley scattering. If the intervalley scattering time is 
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shorter than the relevant time scale for the level spacing 
(i.e. the Heisenberg time), time- reversal symmetry will 
be restored. To probe this idea further, we consider an- 
other observable with a very different time scale, namely 
the conductance of an open cavity for which the time 
scale is the escape time. 

Quantum transport: weak localization — The conduc- 
tance of a cavity attached to two leads [Fig.[TJc),(d)] 
is proportional to the quantum mechanical transmission 
probability from one lead to the other. We use a recur- 
sive Green function method [l6| to find the transmission 
for tight-binding cavities with either abrupt or smooth 
boundaries. First, we focus on the average transmis- 
sion (T ) , where the average is performed with respect 
to the Fermi energy - see Figs.[3lja) and (b) for an exam- 
ple of T{Ef) and its average. Wc find that (T) « 0.5M 
to leading order in the number of open channels in the 
leads M [Fig.[3jb)]: particles are transmitted or reflected 
with about equal probability. The next order correction, 
known as the weak localization correction, is the y-axis 
intercept in Fig.[3jb). It has been studied theoretically 
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for diffusive graphene systems [1^ uM ■ As expected for 
abrupt termination, there is no offset for large enough 
magnetic fields. 

We now focus on the average magnetoconductance of 
our graphene billiards to study the weak localization cor- 
rection in more detail. We compare the magnetoconduc- 
tance data to the semiclassical Lorentzian prediction [2^ : 

(AT(B) } = {T{B) - T(0)) = n/[l + ($o/2^oS)2] (4) 

where TZ is the total magnitude of the effect, is the 
typical area enclosed by classical paths. According to 
random matrix theory (RMT) [ill, IS , 7^ = M/(4M-H2) 
is the difference between the average conductance in sys- 
tems with unitary and orthogonal symmetry (weak local- 
ization is suppressed for unitary symmetry), in agreement 
with the semiclassical theory [23j for large M. 

The numerical results obtained by averaging over 
an energy window are in good agreement with Eq. (jj]) 
[Fig.[3Kc)]. The fit parameter is of the order of the 
billiard area so that weak localization is suppressed 
for a magnetic flux of about <I>o. For the abruptly ter- 
minated billiard with armchair leads, we find 7^ = 0.19 
while the corresponding RMT value is 7?.rmt = 0.20. For 
zigzag leads (in the multi-mode regime), we find 7^ = 0.18 
while 7?,RMT = 0.22. Thus, for the abruptly terminated 
billiards, our numerical results agree with RMT for the 
expected symmetry classes. 

For smooth mass confinement, the expected symmetry 
classes are unitary, both in the absence and presence of 
a magnetic field. Thus, no weak localization correction 
is expected. Numerically, a very small weak localization 
correction is visible: TZ ~ 0.057. We assign the slight 
increase of (AT ) to weak residual intervalley scattering. 

Conductance fluctuations — To show the change in 
symmetry class upon applying a magnetic field for 
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FIG. 3: (Color online) Average conductance: weak localiza- 
tion, (a) Transmission as a function of energy for an abruptly 
terminated billiard [Fig.[ljc)] with zigzag leads (solid line). 
The dashed line shows the number of open channels in the 
leads, M. (b) Average transmission as a function of M for the 
same system with $ = (solid black line, open triangles) and 
$ = 1.6 $0 (dashed red line, full triangles), (c) Change in the 
average transmission as a function of the magnetic flux. Cir- 
cles: Abrupt termination with armchair leads (1-7 open chan- 
nels). The flt (solid line) yields = l.SAs and 7^ = 0.19. 
Triangles: Abrupt termination with zigzag leads (3-7 open 
channels). The flt (dashed line) yields Ao = l.OAs and 
TL — 0.18. Diamonds: Smooth mass confinement [Fig.[ljd)] 
(2-8 open channels). The flt (dotted fine) yields Ao = 0.54As 
and TZ = 0.057 (Parameters of the billiards given in [13]). 



smooth mass confinement, we turn to conductance fluc- 
tuations. Universal conductance fluctuations for the or- 
thogonal symmetry class were found in transport calcu- 
lations on weakly-disordered, rectangular graphene sam- 
ples with zigzag edges [31 . Here, to obtain direct infor- 
mation about the symmetry classes, we investigate the 
magnitude of the conductance fluctuations in chaotic cav- 
ities as a function of energy. The RMT results for the 
variance of the conductance as a function of M are given 
in [2H [Eq. 3(b)] and f22| (Eq. 11), for the cases of the cir- 
cular orthogonal (COE) and the circular unitary (CUE) 
ensemble. 

In Fig. we present the numerical results for the con- 
ductance fluctuations. For the cavities with abruptly ter- 
minated edges, var(r) clearly agrees with the COE result 
when B = 0, while it follows the CUE curve if a magnetic 
field is present. This is as expected from the symmetry 
considerations and weak localization results. 

For smooth mass confinement, Fig.[Hb) shows that the 
magnitude of the fluctuations at zero magnetic field is 
much larger than the COE or CUE values. Rather, it 
is approximately four times the CUE value. When a 
magnetic field is applied, var(T) becomes smaller, about 
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FIG. 4: (Color online) Conductance fluctuations: Variance of 
the transmission as a function of the number of open channels 
in the leads (same cavities as in Fig.[3]l. (a) Abruptly termi- 
nated boundary. B = (black open symbols) and B ^ (red 
full symbols) results are in good agreement with the corre- 
sponding RMT values, orthogonal (COE, black solid line) and 
unitary (CUE, red dashed line). The unitary data uses sev- 
eral values for the magnetic held in the range $ G [0.8, 2.4]"I>o; 
both armchair leads (triangles) and zigzag leads (circles) are 
used, (b) Smooth mass confinement. Zero field (black open 
symbols) and ^ = 2.0 $0 (red full symbols) results are com- 
pared to 1, 2, and 4 times the CUE values (black dotted, red 
dashed, and black solid lines). 



twice the CUE value. This is consistent with the symme- 
try considerations given at the beginning of this paper: 
An ensemble of transmission matrices each with two iden- 
tical blocks implies that var(r) will be 4 times the value 
for a single block. However, an ensemble of transmission 
matrices, each with two uncorrelated blocks, yields the 
sum of the single blocks' values. Since the blocks are 
expected to be unitary in the case of smooth mass con- 
finement, with or without a magnetic field, the result in 
Fig.lHb) follows. 

To summarize, dots formed by mass confinement do 
not follow expectations derived from the effective Dirac 
equation. While the transmission statistics follow from 
the expected block unitary structure, the spectral statis- 
tics show orthogonal or even Poisson statistics. Thus, 
the spectral and transmission statistics follow different 
ensembles! This paradox arises from residual intervalley 
scattering in our system - though the confinement used, 
m(x, y), varies on a scale of 10-30 lattice constants for our 
dots, some weak lattice effects always remain. The time 
scale appropriate for transmission statistics is the escape 
time from the cavity while the time scale for spectral 
statistics is the much longer inverse level spacing. Hence 
if the intervalley scattering time lies between the two, 
different behavior can result. Our study suggests that 
it will be more fruitful to look for smooth confinement 
effects, such as the Berry and Mondragon breaking of or- 
thogonal symmetry without a magnetic field Q , in open 
rather than closed dots. 
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preprint on spectral statistics in nanotube-like structures, 
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